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A CRITICAL ANALYSIS OF THE GRAD APPROXIMATION FOR CLOSING OUT 
THE MAGNETOHYDRODYNAMIC EQUATIONS FOR PLASMAS 


By Willard E. Meador 
Langley Research Center 

SUMMARY 

A critical analysis is made of the utility of Grad's 13-moment distribution functions 
in closing out the generalized magnetohydrodynamic equations for plasmas. In particu- 
lar, the inaccurate particle velocity dependence of the Grad approximation is shown to 
yield significant errors in the collisional transfers of such quantities as momentum and 
energy as functions of the moments (variables of the problem) considered. These errors 
directly affect the relations between moments which are derived by means of the afore- 
mentioned closing-out process; consequently, when the Grad functions are written so as 
to incorporate the moment relations, considerable doubt is cast upon the ability of such 
functions not only to describe the original 13 moments but also to predict accurately other 
plasma properties (e.g., the entropy) as well. Numerous comparisons with the exact 
first- and second-order perturbation solutions of the Boltzmann equation are made in 
order to substantiate these conclusions. A study also is made of the convergence of 
Grad-type functions as more moments are added, in consequence of which at least 
16 moments are often necessary. 


INTRODUCTION 

Macroscopic equations of change for such velocity moments of the Boltzmann 
kinetic equations as constituent number densities, mean velocities, temperatures, 
stresses, and heat fluxes in plasmas have been generated by many authors, but with 
widely varying accuracies with regard to the net effects of interparticle collisions. A 
particular example of the more rigorous approaches is the work of Everett (ref. 1), in 
which Grad's 13-moment velocity distribution functions are utilized for closing out the 
moment equations, that is, for expressing the collisional transfer terms as functions of 
the moments (variables) of the problem. When coupled with Maxwell's electromagnetic 
field equations, the resulting magnetohydrodynamic equations constitute a closed set. 

The efficacy of such developments obviously depends on the accuracy of the distri- 
bution functions employed in the closing-out process. In particular, errors in the 
detailed dependence of the distribution functions upon particle velocities may lead to 


serious errors in the relations between the collisional transfer terms and the moments 
of the problem and also in the relations between moments which are deducible from the 
final closed set of macroscopic equations. Examples of moment relations include the 
specifications of heat fluxes and pressure tensor elements in terms of the electron diffu- 
sion velocity. 

A prime purpose of the present research is to investigate the accuracy of the 
moment relations obtained by using the Grad 13-moment velocity distribution functions in 
the aforementioned closing-out process. This is accomplished by first assuming the dis- 
tribution functions to be expressible as perturbation series, then employing the closed set 
of macroscopic equations to compute the heat flux, a certain higher moment, and the 
second-order contributions to the pressure tensor as functions of the electron diffusion 
velocity, and finally comparing the results with those obtained from exact first- and 
second-order perturbation solutions of the Boltzmann kinetic equations. The character 
of each perturbation series is such that the zeroth-order term corresponds to a local 
Maxwellian distribution function, whereas the first- and second-order terms correspond 
to corrections which are linear and quadratic, respectively, in the electron diffusion 
velocity or in its driving force. 

Although this form of solution is not required by the Grad N-moment method, which 
is a more general approach that includes the perturbation expansion as a special case, 
there are two overriding advantages of the present procedure. They are, first, that such 
solutions are physically meaningful and, second, that exact solutions of the corresponding 
integro-differential perturbation equations are readily found if modifications to Meador f s 
collision model (ref. 2) are substituted for the Boltzmann collision terms. Thus, all ref- 
erences herein to the order of a function or quantity indicate the highest power retained 
of the electron diffusion velocity or its driving force. 

A study also is made of the convergence properties of Grad N-moment functions, 
whereupon it is observed that at least 16 moments are often necessary for the closing-out 
method to yield accurate moment relations and that the addition of higher first-order 
moments to the Grad functions is equivalent in a systematic fashion to the calculation of 
higher Sonine approximations to the Chapman-Enskog distribution functions (ref. 3). 
Several types of interparticle interaction potentials are considered, including those of the 
real (as opposed to Lorentz) fully ionized gas. 

Another objective of the present research is to investigate the ability of Grad 
N-moment functions to predict accurately additional plasma properties (e.g., the entropy). 
This particular study involves the substitutions of the following two sets of moment rela- 
tions into the original Grad functions: (1) the exact moment relations obtained from solu- 

tions of the Boltzmann kinetic equations and (2) those obtained from the closing-out pro- 
cess. The accuracies of the entropies thus calculated are found to be highly sensitive to 
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the plasma conditions; more specifically, the errors are especially large when the elec- 
tric field and the temperature gradient are so related as to yield zero heat flux with 
respect to the average electron motion. 

Finally, the Grad 13-moment functions are modified through second order by first 
imposing entropy maximization, which is the approach of information theory, and then 
referring the electron particle velocity to the electron frame of reference, which is the 
method suggested by Everett (ref. 1). Although both modifications lead to improvements 
in the evaluation of the pressure tensor, the one proposed by Everett gives the best 
agreement with the exact calculations. 

Unless otherwise noted, the following simplifying assumptions apply throughout the 
present report: time independence for all parameters, infinite mass for the heavy parti- 

cles which are considered fixed scattering centers for the electrons, zero plasma flow 
velocity, zero applied magnetic field, and spatially homogeneous fluxes, temperatures, 
pressures, and all other quantities through the second perturbation order. 

SYMBOLS 

A(£) ratio of impact parameter integrals 

b impact parameter 

B magnetic field 

B unit vector in direction of magnetic field 

c e electron particle velocity 

cfj heavy-particle particle velocity 

e magnitude of electron charge 

E electric field 

f e electron distribution function 

f^ Maxwellian contribution to electron distribution function 

fj heavy-particle distribution function 


3 



hjo 


gi(y)>g2 (?) 

AAA 

T 

k 


P e 

p e,xy 

o 

o 

P 

*e,xy 
e,xy 

Qi 

Q2 

s 

s<°) 


functions defined in equations (81), (84), and (85) 
indices 

unit vector in x-, y-, and z-direction, respectively 

electron current density 

Boltzmann's constant 

mass of electron 

electron number density 

number density of heavy particles 

electron partial pressure 

xy-component of electron pressure tensor (subscript used indicates appro- 
priate component) 

traceless electron pressure tensor 

xy-component of traceless electron pressure tensor referred to plasma 
motion (subscript used indicates appropriate component) 

xy-component of traceless electron pressure tensor referred to electron dif- 
fusion velocity (subscript used indicates appropriate component) 

heat flux vector, £1 - -f /3n (see table I) 

vector, ^2 " + - 4 - Po ( see table II) 

integral defined by equation (16) 

entropy density 

equilibrium entropy density 
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Sj 

t 

T e 

n 

U 

v e 

x 

x,y,z 


collisional transfer of reduced flux j3j 
time 

electron temperature referred to plasma motion 
electron temperature referred to electron diffusion velocity 
unit tensor 

electron diffusion velocity 
integration variable 
Cartesian coordinates 


0 fe 
at 


I time rate of change of 


at 

' 8 if 

at 



contribution to 

f 8f e 

C,<Pl 


\8ty 

1 

c,<?2 

contribution to 

(*e 

\ at/ 


f e due to collisions 


from first-order perturbation function 


from second-order perturbation function 


“1 >“2 

y 

6 

e 

V 

e 


constants 

velocity moment defined by equation (3); for example, /3 q is reduced elec- 
tron diffusion velocity and j3j is reduced heat flux 

reduced electron particle velocity 

variational operator 

azimuthal angle for electron — heavy-particle collision 
heat flux parameter defined by equation (94) 
polar angle in spherical coordinates 
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fij Lagrange multiplier 

£ effective interparticle interaction parameter 

a electrical conductivity 

T g collision time for entropy production 

r a collision time associated with electron diffusion 

cp azimuthal angle in spherical coordinates 

<Pl first-order perturbation function 

cp 2 second-order perturbation function 

cp e electron perturbation function 

X scattering or deflection angle in electron — heavy-particle collision 

co cyclotron frequency 

Subscripts: 

e electrons 

j heavy particles 

x,y,z x- ? y-, and z-components 


Special notations: 

<> average over velocity space 

<>n average over velocity space using only Maxwellian distribution 




factorial employing only odd numbers; for example, (5).'! = 5- 3- 1 = 15 


Primed quantities in collision integrals denote conditions after collisions, as 
opposed to unprimed quantities which denote conditions before collisions. 
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MACROSCOPIC EQUATIONS OF CHANGE 


A convenient initial step in the investigation of first-order moment relations is the 
derivation from the macroscopic equations of change of a recursion-type formulation of 
the collisional transfer terms. Such is the purpose of the present section, in which the 
fundamental working equations are deduced. Second-order moment relations are dis- 
cussed in a subsequent section. As explained in the Introduction, the restriction to first 
order is imposed by neglecting squares and higher powers of the electron diffusion veloc- 
ity (or its driving force). 

The macroscopic equations of change for electrons are generated by taking veloc- 
ity moments of the Boltzmann integro-differential equation (ref. 3) 


Me f 2k T e \ 1/2 - 

ir + v^w y 



( 1 ) 


in which (8f e /8t) c represents the set of collision integrals and y is the reduced elec- 
tron particle velocity defined by 


y = 



( 2 ) 


There is one such equation of change for each of the moments regarded as variables of 
the problem. 

As an example of this procedure, the equation of change for the moment 


/3j s (y^y) = n e -1 J y 2 jy i e dc e (3) 

is obtained from equation (1) by multiplying each member of that expression by c^(T e , 
integrating over velocity space, and employing vector calculus and Gauss 1 divergence 
theorem in the usual manner (ref. 3). The result is 

+ ^(^rJ(< r2j ) S + 2jS • <r 2j - 2 rT)) 

(4) 
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where Sj refers to the collisional transfer of /3j and is expressed as 



( 5 ) 


It is evident from equation (4) that if 0, 1, and 2 are the only values considered for 
the index j, the variables of the system (in an unreduced sense) will be comprised of the 
28 moments corresponding to n e , /3q, /3j, and the symmetric tensors ^yy^ , ^y2yy^, 

and ^/4yy^ . Likewise, j selected as 0 and 1 involves the 19 moments corresponding to 
n e , /3i, and the symmetric tensors ^yy*^ and <^y2yy^, which set can be reduced to the 

familiar 13 moments of Grad (ref. 1) by using a distribution function to express <y2yy*^ in 
terms of the remaining variables. Such effort is part of the closing-out process, as is 
the elimination in a similar manner of the collisional transfer term Sj in equation (4). 
The number of moments ultimately retained as variables of the system is obviously a 
matter of some choice, there being no strict limitations on where the closing -out process 
must stop. 


Because of the previous restrictions to spatial homogeneity and time independence, 
the first two terms on the left side of equation (4) are zero. In addition, since the elec- 
tric current is assumed to generate the only magnetic field and the product /3j x B is 
therefore second order by reason of Maxwell's equations, the first-order form of equa- 
tion (4) can be written 



The zero subscripts refer to velocity averages taken with respect to the zeroth-order 
Maxwellian distribution function given by 

m e \ 3 / 2 _y2 


f e n e 


27rkT ( 


(7) 


Equation (6) is further reduced by utilizing 


<Jy 2 i) 0 = ng 1 j* y 2 34 0 W e = 2“j(2j + 1).’.’ 


( 8 ) 
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to write the equation of motion 


and the recursion relation 

S j+1 = |(5 + 2j)Sj (10) 

Equation (10) depends only on the assumptions of first-order theory and spatial and time 
independence. This equation is the working equation for the computation of first-order 
moment relations, together with equation (5) and a formulation of (9f e /3t) c yet to be 
given. 

More specifically, if an electron distribution function of the Grad N-moment variety 
is employed in equations (5) and (10), all moments /3j can be expressed in terms of the 
reduced electron diffusion velocity /3q. Hence, /3g is regarded as input data and thus 
is independent of the approximation to the distribution function. Equation (9) is of use 
only when the dependence of j3j upon the driving electric field is desired. 

Grad N-Moment Distribution Functions 

A convenient way of writing Grad's 13-moment electron distribution function for the 
purpose of closing out equations (5) and (10) is as follows (ref. 1): 

f e = f e 0) (l + <?e) (ID 

where 

cp e = 2P 0 • y + - |j(0i y + l-P e :yy (12) 

0 

and P e is the traceless pressure tensor defined by 

Pe = 2P e (<^ - I ®) " TQT $ 4 0> We d? e (13) 

Although equation (12) and its first-order form (see appendix A for an alternate 
derivation) 

<Pe = 2 % ' y + \ (r 2 ~ f)(/3l -!&)•? (14) 



( 9 ) 
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produce identities for the moments when the proper integrations are performed, this 
characteristic is not sufficient to guarantee accurate values of Sj when equations (11) 
and (14) are employed in equation (5). The integrals defining these moments are quite 
different from those defining the Sj, so that errors in the detailed y-dependence of equa- 
tion (14) will be reflected in the results from equation (5) and also in the moment rela- 
tions subsequently derived from equation (10). 

Before any distribution function can be employed in equation (5), however, the colli - 
sional derivative (8f e /9t) c must be expressed in terms of <p e . This is accomplished 
by using the following simplified first-order collision model developed by Meador (ref. 2): 

= v l-(V{),(o) (15) 

c m e aR 04 



in which Rjj is the integral 



and | is an effective interaction parameter defined in such a way as to guarantee the 
correct relation between the heat flux and the electron diffusion velocity on the basis of 
the exact solution to the corresponding Boltzmann equation. 

Reference 2 further shows that equation (15) is entirely adequate for describing a 
great variety of realistic plasmas, even though to the first Chapman-Enskog perturbation 
order the description is exact only for a Lorentz gas (i.e., negligible electron-electron 
interaction effects either in a slightly ionized gas or in a fully ionized gas having large 
ionic charges). The parameter £ in the latter cases is the exponent in the electron — 
heavy-particle inverse-power interaction potential. 

The substitution of equations (14) and (15) into equation (5) yields 

1 - 2j)£ + 

which relation combines with equation (10) and the identity 



Sj( 13 -moment) = .! Vl3^| 

2m„CTR^.£ ^ 


2m e aRQ^t 


to give 


R -l,2j+7 " ^ -1 |] 3 + ^ " ?l R -l,2j+5 
^(13-moment) 

2[(13 + 6j)£ 2 - 8(2 + j)4 + 16] 


(18) 

( 19 ) 
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In particular, when j = 0, the expression 


/S* (13 -moment) = ^(ng 2 - ,g| + l g) ^ ( 20 ) 

2(l3£ 2 - 16£ + 16) 

represents the first-order culmination of the present method for closing out the macro- 
scopic equations of change at the 13-moment level. 

Likewise, for the other Grad functions of present interest, the substitution of 

<p e (5-moment) = 2/3 q • y (21) 

and equation (15) into equation (5) gives 

cT/ c .v e2n ! R 13 R -l,2j+5 -? (00 , 

Sj(5-moment) = J — /3 q (22) 

m e aR 04 


and the substitution of equation (A5) and equation (15) into equation (5) gives 


Sj(16-moment) = - 


e 2n l R 13 R - 


l,2j+5 


2m e aR 2 4 ^2 


i 


(3 - 8j + 4j 2 )^ 2 + 16(1 - j)£ + 16 


00 



+ 4j - 4j2)^2 . 8(1 _ 2 j)| - 16 


01 



4j 2 )^ 2 + 16j| 



Accordingly, equations (10), (18), and (22) yield 


(23) 


5 + 2j _ R -l,2j+7 _ (3 + j)| - 2 
2 R -l,2j+5 £ 

so that 


(24) 


|(5-moment) = 4 (25) 

In a similar manner, the following expression is deduced from equations (10), (18), 
and (23): 
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/32(16-moment) = 35 


(23 + 56 j + 20j 2 )? 3 - 4(27 + 20j + 4j 2 )? 2 
+ 16(5 + 4j)? - 64 ] ^ (21 - 32j - 20j 2 )? 3 

+ 4(l9 + 12j + 4j2)^2 _ i 6 (i + 4j)| + 64j^ 0 
- - 44 j - 20j 2 )^ 3 + 4(l9 + 16] + 4j2)^2 


- 16(3 + 4j)? + 64 


/3j(16-momenty 


(26) 


Equation (26) is further simplified by solving simultaneously the j = 0 and j = 1 
expressions, which are 

1 


/32(16-moment) = 


5^21? 3 + 76? 2 - 16? + 64)is 0 - 4(3? 3 + 76? 2 - 48? + 64)0 1 (16-moment)J 


4(23| 3 - 108£ 2 + 80£ - 64 


(27) 


and 


g 2 (16-mon,ent) - ’fc 61 * 8 - 156{2 * 112{ - ««)%(16-mo m e..) - sfrljf - »of j 80 t ' 64 fe] ( 28 ) 

4 ( 99| 3 - 204| 2 + 144 | - 64 ) 

to obtain the first-order result 


/3i(16-moment) 


5(349? 4 - 416? 3 + 672? 2 - 512? + 256) j 
2(425? 4 - 816?3 + 1120? 2 _ 768£ + 256) 


(29) 


Only one conclusion can be drawn from the differences between equations (20) 
and (29): In spite of the fact that the macroscopic variables /3q and /3i explicitly 
appear in each of equations (14) and (A5), the 13 -moment and 16 -moment distribution 
functions do not predict the same expressions for the collisional transfer terms 
Sj^/3o,/3i). Hence the accuracy of the closing-out procedure may be quite sensitive to the 
number of moments considered as variables in the original cp e . The one exception 
occurs for Maxwellian molecules (? = 4), in which case all results agree with 
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(30) 


^(5-moment) = n^ 1 J y^e’^O • y)dc" e = | 0 O 

and 

/3 2 (5-moment) = n^ 1 j* y 4 ^ o ^2^ 0 ’ y)dc e = ^ @0 = 7 (Pl ~ f Poj ( 31 ) 

because of the relation 

(p e (13-moment,£=4) = <p e (16 -moment, !=4) = cp e (5 -moment) (32) 

In analogy with equation (25), the combination of equation (20) with the j = 1 form 

5(l7£ 2 - 16| + 16) - 

0i( 13 -moment) = -j y /3 0 ( 33 ) 

2(l9! 2 - 24! + 16) 

of equation (19) gives 

^(13-moment) = 4 or 4/3 (34) 

The bounds of the closing-out procedure are exceeded in both instances by the application 
of equation (10); in particular, the macroscopic equation of change for 0j with j = 1 is 
used in deriving equation (25) and the macroscopic equation of change for 02 is used in 
deriving equation (33), whereas the accurate descriptions of the 5-moment and 13-moment 
distribution functions stop (unless £ = 4) with 0 q and 0j, respectively. A definite 
upper limit on the index j is always imposed in equation (10) if the collisional transfer 
terms are computed from a Grad N-moment function. 

Comparisons of Heat Fluxes 

Reference 2 shows that the exact (or infinite moment) first-order distribution func- 
tion corresponding to an applied electric field and the collision model of equation (15) can 
be written 


^ e (°°-moment) 


2R 04 

r 13 


y (4/£)-lj3 0 . y 


Accordingly, 


^(oo-moment) = 2!"1(! + 1)/3 q 


(35) 


(36) 
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and 


^("-moment) = 2£~^3| 2 + 5| + 2)j3g (37) 

from equations (3) and (35). 

The following summary of heat fluxes obtained from equations (30), (20), (29), 
and (36) can thus be used to assess the convergence properties of Grad-like distribution 
functions in closing out the macroscopic equations of change: 


/3j(5-moment) - (3q = 0 

Li 

^(13 -moment) - | js 0 = - — i3 0 

2 134 2 - 164 + 16 

3 ,+« ^ 104(4 - 4)(l94 2 - 244 + 16) - 

/3 1 (16-moment) - x % = 7 ; L 00 

4254 4 - 8164 3 + 11204 2 - 7684 + 256 

and 


(38) 

(39) 

(40) 


/3j( 00 -moment) - | P Q = p Q (41) 

Numerical calculations based on equations (38) to (41) are given in table I for effec- 
tive interparticle interaction potentials corresponding to a Lorentz fully ionized gas 
(£ = 1), a real fully ionized gas including electron-electron collisions (£ = 1.6674 from 
ref. 2), a Lorentz gas with £ = 2, a Lorentz gas composed of Maxwellian molecules 
(£ = 4), and a Lorentz gas of rigid spheres (| = °°). 

TABLE I.- NONDIMENSIONAL HEAT FLUX Qj = 01 - - iSq 

2 


N 


Ql I&0 for 4 of 

- 


1 

1.6674 

2 

4 

OO 

5 

0.000 

0.000 

0.000 

0.000 

0.000 

13 

1.154 

.764 

.556 

.000 

-.385 

16 

1.521 

.704 

.507 

.000 

-.447 

OO 

1.500 

.699 

.500 

.000 

-.500 


Although the errors in the heat flux at the 13-moment level are not especially serious in 
first-order calculations, they may well be large compared with the higher order 


14 




corrections (to linear flux theory) discussed by Everett. In addition, the errors in 0j 
are fed back to the distribution functions themselves and will affect the strongest 
y-dependent terms therein; consequently, the errors may be compounded in other types 
of velocity integration (including the macroscopic equations of change for higher order 
quantities). 

Another demonstration of the fact that the Grad function is less accurate in the 
present context than was commonly supposed is obtained from similar comparisons for 
02- The expressions analogous to equations (38) to (41) are quite complicated, however, 
especially for the 16-moment velocity distribution; thus, it is more convenient in that 
case to compute the fluxes directly from equations (28) and (29). There results 


02(16-moment,|=l) - 7|0j(16-moment,^=l) - 0 q| = 0.96800 

02(16-moment, £=1.6674) - 7|0i(16-moment, £=1.6674) - ~ 0qJ = -O.2870Q 
02(16-moment,|=2) - 7|0j(16-moment,|=2) - 0qJ = -O.3230Q 
02(16-moment,£=4) - 7 |0i(16-moment,(;=4) - 0 qJ =0 


(42) 

(43) 

(44) 

(45) 


and 


02(16-moment,£=°°) - 7 


0l(16-moment,|=°°) 


| 0q] = 0.49400 


(46) 


The remaining examples are not as difficult to handle and the following equations 
may be derived from equations (3), (14), and (39) for the 13-moment function and from 
equations (31), (36), and (37) for the 5- and ^-moment functions: 


02(5-moment) - 7 0i(5-moment) 


& 


- 7 (^: 


j32(13-moment) - 7 0^13-moment) 


•!?o]=o 


and 


/3 2 (°°-moment) - 7 ^(^-moment) 


I? 


£1 . (I - 4)(3« -4) 


-| ? o " 




(47) 

(48) 

(49) 
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Numerical calculations based on equations (42) to (49) are given in table II, where 
it is seen that the 13 -moment analysis is completely inadequate (except for Maxwellian 
molecules) and even the 16 -moment analysis yields significant errors when compared 
with the exact answers. 


TABLE II.- NONDIMENSIONAL FLUX Q 2 = j3 2 - 7^ + j3 0 



Entropy Calculations 

Another important application of the present theory involves calculations of the 
entropy density difference (ref. 2) 

s (o) - s =|J4°Vi d c e (50) 

This and related quantities provide excellent criteria on the utility of the closing-out pro- 
cedure for predicting plasma properties other than the basic moments appearing in the 
Grad-like functions. 

It is further shown in reference 2 that if the collision time r s for entropy produc- 
tion is introduced according to the definition 


_ t g ( s (°) - s) 

l — - - ■— 


n e k ^o 


(51) 


where t q is the familiar collision time 


e*n P m^cr 


(52) 


the generalized Ohm T s law can be written (through linear terms in cor s ) in a form com- 
pletely analogous to the mean-free-path result. Hence, 


j = <?(i 


E + WT S B x E 


(53) 
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The computation of r s from equations (50) and (51) is thus equivalent to finding the Hall 
conductivity in the case of small magnetic fields. 

The evaluation of equation (50) requires the following set of perturbation functions 
obtained from combinations of equations (21), (14), (A5), (35), and (38) to (46): 


cp e { 5-moment) = 2/3 q • y (54) 

. „ 4(|(4 - £)y2 + 9^ 2 - 18£ + £1 - _ . BB . 

<p p (13-moment) = — — /3 n • y (55) 

13| 2 - 16 £ + 16 

<p e (16 - moment, £=1) = JL(o.968y 4 + 3.871y2 _ 0.645)/3 0 • 7 (56) 

cp e (16-moment, £=1.6674) = -JL(o.287y 4 - 6.937y 2 - 2.669^ 0 ‘ 7 (57) 

<p e (16-moment,£=2) = - — (o.323y 4 - 5.807y 2 - 5.807)/3 0 • y (58) 

35 

<P e (16-moment,i;=4) = 2/3 q • y (59) 


<p e (16-moment,|=°°) = — fo.494y4 - 6.588y^ + 29.65^/3 q • y 

35 v • 


and 


cp e ( 00 -moment) 


2R 04 (4/|)-lx - 

- 5 — y &o ' y 

R 13 


(60) 


(61) 


Numerical calculations based on equations (50), (51), and (54) to (61) appear in 
table HI for a variety of effective interparticle interaction potentials and number N of 
moments in the Grad-like distribution. Except for the trivial case of £ = 4, the 
13-moment analysis yields an error in ^ T s/ T a) " ^ that ran S es from 17 percent for 
£ = 1.6674 to 67 percent for £ = 
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TABLE HI.- RATIO OF COLLISION TIMES r g AND T a 


N 


T s 

/ t ct for £ of 

— 


1 

1.6674 

2 

4 

oo 

5 

1.000 

1.000 

1.000 

1.000 

1.000 

13 

1.533 

1.233 

1.123 

1.000 

1.059 

16 

1.979 

1.203 

1.109 

1.000 

1.094 

CO 

1.933 

1.199 

1.104 

1.000 

1.178 


The values of T s / T 0 for £ = 1 coincide exactly with the first, second, third, and 
infinite Sonine approximations (ref. 2), respectively, as found from solutions of the 
Boltzmann equation. That this is no accident follows from the observation that equa- 
tion (55), for example, is identical for £ = 1 with the second Sonine approximation 
derived in appendix C of reference 2. It is not immediately obvious why Grad-like func- 
tions reduce to Sonine approximations to the solution of the integro-differential kinetic 
equation when the higher moments /3j are related to /3 q through the closing out of the 
macroscopic equations of change; however, a key factor may be the appearance of Sonine 
polynomials (in y) in equation (A5). In any event, the fact that the addition of higher first- 
order moments to the Grad function in the present method corresponds in one-to-one fash- 
ion to higher Sonine approximations to the Chapman-Enskog distribution function is of 
considerable interest because it provides an alternate and equivalent method of solution. 

One final illustration of this equivalence is given in table IV, where the ratio of t s 
to t g computed in the present paper is compared for the real fully ionized gas 
(| = 1.6674) with the Sonine approximations (denoted by N') of reference 2. The small 
differences that do prevail at each level are caused by the use of the collision model of 
equation (15) in the one case and the more rigorous Chapman-Enskog treatment of 
electron-electron encounters in the other. Oddly enough, the much simpler approach 
yields slightly better values with respect to the exact “-moment numbers. 

TABLE IV.- RATIO OF COLLISION TIMES T g AND t q FOR 
A REAL FULLY IONIZED GAS (£ = 1.6674) 


N 

T s/ T a 

(Grad-like) 

N’ 

T s/ T a 

(Chapman-Enskog) 

5 

1.000 

1 

1.000 

13 

1.233 

2 

1.259 

16 

1.203 

3 

1.207 

CO 

1.199 

CO 

— 
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Other Applications of Grad-Like Functions 

It has now been shown that Everett's method of utilizing the Grad 13-moment dis- 
tribution function for closing out the macroscopic equations of change leads to substantial 
errors in the relations between /3 q and the higher fluxes These errors, in turn, 

are fed back to the distribution function itself, so that subsequent calculations of quantities 
such as the entropy are similarly affected. What is not clear, however, is whether the 
errors in the entropy are caused primarily by the approximate forms of the distribution 
function or by the use of incorrect moments therein. Even though such distinctions are 
improper in the strictest sense, the possibility that Grad-like functions employing exact 
/3j -relations will yield improved entropy values should be investigated. 

Two cases are considered for this purpose: (1) the spatially homogeneous Lorentz 

fully ionized gas (£ = 1) already discussed; (2) a fully ionized Lorentz gas in which there 
is no heat flux relative to the average electron motion and which implies a temperature 
gradient in fixed relation to the electric field. Only in case (1) does the 13-moment 
function prove to be reliable. 

Since equations (41) and (49) yield 

%-|% = |?0 <« 2 > 

and 

h ' 7 (?i - f %) - 1 % < 63 > 

for case (1), the Grad-like perturbation functions may be written as follows from equa- 
tions (21), (14), (A5), and (35): 


<p e (5-moment) = 2/3q • y 

(64) 

<p e (13-moment) = y 2 - ljjsf 0 • y 

(65) 

<P e (16 -moment) = ^ y 4 + | y 2 - ^Jp Q • y 

(66) 

1/2 

<p e (°°-moment) = ?-j— y 3 /3 Q • y 

(67) 
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Numerical calculations based on equations (50), (51), and (64) to (67) are presented 
in case (1) of table V, where it is evident that a major improvement is obtained over the 
corresponding 13-moment result in table m. What small error remains is thus attribut- 
able entirely to the form of the distribution function rather than to an incorrect relation 
between [3^ and /3q. 

TABLE V.- RATIO OF COLLISION TIMES r s AND r CT FOR 
A LORENTZ FULLY IONIZED GAS 


Conditions 

N 

7 s/ T o 

Case (1): 

5 

1.000 

No gradients, 

13 

1.900 

h = 4 ?0 

16 

1.932 


oo 

1.933 

Case (2): 

5 

1.000 

Temperature gradient, 

13 

1.000 

% = 5 %/ 2 

16 

1.804 


OO 

1.865 


Case (2) requires the solution of equation (1) with the temperature gradient included, 
which is available from reference 2 in the form 


1/2 

» (-moment) =- I y (V9-l 

e Rl3 V m e / 


ar i+ h 2 -!i vl " Te 


y (68) 


for zero pressure gradient. This expression can be used to evaluate |3^ and /3 q as 


/3^ = 7T“3/2 ^ e~r-)2y(p e d y = - 


U+ l) To/ nl e \ 1/2 ro„^ 


(r m e 


2kT c 


2e|E + k(^+ 4)VT e J (69) 


and 


^0 _ 





(70) 
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the combination of which with equation (68) and 




( 71 ) 


finally gives 


(P e (°°-moment) = 


r 3 (3y 2 - 20)/3 q • y 


if k = 1. 

Since the appropriate moments of equation (72) yield 


% ' 7 (^1 ' 4 ^o) " “ T ^0 


(72) 


(73) 


the following set of perturbation functions are obtained when equations (71) and (73) are 
imposed on equations (21), (14), and (A5): 


c? e (5-moment) = cp Q (13 -moment) = 2/3q • y (74) 

and 

c> e (16-moment) = - |^y 4 - 7 y 2 + * Y (75) 

Numerical calculations based on equations (50), (51), (72), (74), and (75) are given 
in case (2) of table V, where it is seen that the 13-moment approximation to the ratio 
of Tg to t q is inadequate. The reason is obvious: Since cp e (13 -moment) is the 
same as <p Q (5 -moment), and since f e (5-moment) corresponds to a Maxwellian distri- 
bution referred to the electron frame of reference, the 13 -moment function is unable in 
this case to describe the important frictional effects arising from electron-ion collisions. 
At least the 16 -moment approximation is needed to safeguard against this occurrence. 

Finally, one might consider the prediction of by the 13 -moment distribution 
function. Since the corresponding integration of equation (14) yields 


?2-’(?l-|?o) = ° < 76 > 

regardless of the relation between and /3 q, there can be no agreement at this level 
with equation (63) or (73). Serious difficulties are therefore unavoidable in some appli- 
cations of the Grad function, even when the best values of the fluxes are employed. 
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SECOND-ORDER THEORY 


Exact and Grad Solutions 

As an example of the utility of the closing-out method for predicting higher order 
plasma properties, the second-order traceless pressure tensor so determined is com- 
pared in the present section with calculations involving a direct second-order solution of 
the Boltzmann equation. Errors in the results obtained by using the Grad 13-moment 
approximation are expected to arise from the fact that the pressure tensor term in equa- 
tion (12) incompletely describes higher order contributions in much the same manner that 
the electron diffusion velocity and heat flux incompletely describe first-order effects. 

Since an exact second-order solution is necessary in this investigation, the physical 
problem should be as tractable as possible and simpler than the one considered by 
Everett (ref. 1). The essential requirements are satisfied quite well by a hypothetical 
plasma in which an electric field is applied in the z-direction, the magnetic field is zero 
(in violation of Maxwell's equations), there are no spatial variations of any macroscopic 
quantity, and a Lorentz-like collision model is applicable. 

The pertinent second-order perturbation equation can be written from equation (1) 
as follows: 



% 9 

r a dy z 




(77) 


Primed symbols refer to quantities after a collision (as opposed to unprimed symbols for 
quantities before a collision) and cp e is split into additive and ^ parts, where 

cp g is second order and 


2R 04 (4/£)-l fl 
<p i = — — r /3 0 r z 


R 


13 


(78) 


from equation (35). Also used in the derivation of equation (77) was Ohm's law and 
equation (52). 

It is evident from the energy balance expression (ref. 3) 


3 8 »e = r 2 Vo 

2~ E,J ~rr 


(79) 


that the time derivative of the electron pressure is second order. Moreover, since the 
number density n e is time independent by reason of the electron equation of continuity, 
there results 
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( 80 ) 


!^,(o)L 2 3\ al "Pe .4P4(, 3 

at e V 2/ at 3r a V 2 


from equations (7) and (79). 

A physically meaningful solution is obtained by substituting the trial function 

2 .2\ 


<P 2 = [^z " r ) g l( y ) + 

into equation (77) and using equation (80), the relation 


9 /,(o) .. \ _ 2 4° >R 04' 3 0..(4/e-3 


9y, 1'rnl " 


R 


13 


y 


- 2|y 2 - 4^i)r 2 


and the collision integral 

K* ~ 4 2 ) 


b db de = 


A(£)R 


13 


m £ 


,1/2 


2n j T a R 04\ 2kT e 
derived in appendix B. The resulting equations are 

*2 .2 


r 


-4A 


( 3 ''z - y 2 ) 


M _ 8R 04% r (8/e-4/..2 _ 
2 


9A(?)Rf 3 


y 


2£ 


and 


g 2 (y) 


4 # 

R 04 

3r a 

R 13 


y' 


y‘ 


y + 


Since neither the electron number density 


ne=j’4 0)dE e + j’4° > ?'2 dir e = r >e 


nor the electron temperature 


m e 

T p = — - 
e 3n„k 


(j* c e4° )d? e + J c e4°^2 d? e) = T e 


(81) 


(82) 


(83) 


(84) 


(85) 


(86) 


(87) 
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is affected by the following combination of equations (78), (81), (84), and (85) is the 
correct perturbation function through second order: 


_ 2R 04% (4/8-1 8R 04^ (8/Q-4/ 2 4 - j\L 2 

6 R 1S ^ 9A( S )Rf 3 V 24 )\ ** 



+ 


4$ 

R 04 

2t ct 

R 13 


y (4/0-l 




-r 2 + | 


( 88 ) 


The traceless electron pressure tensor is determined by the substitution of equa- 
tion (88) into equation (13), whereupon one obtains 


and 


o 

P 

e,xy 


o o 

P = P 
e,xz e,yz 


= 0 


o o o o 

P = -P - P = -2P 
e,zz e,xx e,yy e,xx 


64(^ + 1 ) R 04 R 22 „ 0 2 
n PpPn 
45^A(^)r2 3 


(89) 


(90) 


Equations (89) and (90) are the exact second-order tensor elements required for 
comparisons with the approximate Grad 13-moment predictions. It is essential, how- 
ever, that the same physical problem be considered in all calculations; accordingly, the 
application and closing out of the macroscopic equation of change for P must corre- 
spond to the present plasma conditions, which are not the same as those of Everett. 

The multiplication of each term of equation (1) by m_c* and subsequent integra- 

e e , z 

tion over velocity space yields the following P g zz expression 



in which Chapman and Cowling's concepts (ref. 3) of inverse collisions and microscopic 
reversibility, together with several of the assumptions in the last paragraph of the 
Introduction, are employed in the successive simplifications of the collision integral. 
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Further simplifications of equation (91) are obtained by the use of Ohm's law and 
equations (79) and (83). There results 

2 9AU)R irj (■ 2 l-(4/i;)/„ 2 2^ . 

% = 64 T R§. J r V !y z ' 7 H ^ < 92 > 

where the cf> e is that of equation (12) if the unmodified Grad 13-moment function is used 
for closing-out purposes. As mentioned previously, substantial errors may arise from 
the fact that the pressure tensor term in equation (12) incompletely describes second- 
order contributions in much the same manner as the electron diffusion velocity and heat 
flux incompletely describe first-order effects. 

All such second-order corrections to equation (12) are available through straight- 
forward extensions of the information theory outlined in references 4 and 5 and appen- 
dix A. More specifically, if the entropy is maximized subject to constraints on n e , 

__ __ 

Jq, /3p and P e , which are interpreted as the observer's complete first-hand knowledge 
of the plasma, the following modified 13 -moment perturbation function is obtained: 




= 2(Wf rir 2 * l - T)\ + p e,zz[ 3 rl - r 2 ) * 




(93) 


The parameter r] introduced in this expression is defined by 


1 


^0 


(94) 


A difficulty with equation (93) is apparent from the fact that the corresponding 
entropy is associated more with the observer's uncertainty about the plasma than with 
the disorder predicted by the exact solution to the Boltzmann equation. If the number of 
constraints employed in the entropy maximization is too small, the observer's uncer- 
tainty may be so unreasonably large that the two correction terms in equation (93) will 
substantially overcompensate the errors inherent in the unmodified Grad approximation. 
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The substitution of equation (93) into equation (92) finally yields 


e,zz 


(modified 13 -moment) = 


40|Rq4 


9(3? - 2)A(?)R 13 R _ 1j5 


4?7 

75 ? 2 


(7 ? 2 - 32? + 16 ) 77 + 10?(? - 4) 


’ l V 3 0 


(95) 


As an examination of this development demonstrates, all that is necessary for the conver- 
sion of equation (95) to the unmodified Grad 13-moment result is the formal replacement 
of 77 with zero. Calculations of the latter type are presented in the second column of 
table VI for several electron — heavy -par ticie interaction potentials and may be compared 
with the exact results obtained from equation (90) and listed in the fourth column. 


TABLE VI.- TRACELESS PRESSURE TENSOR COMPONENT 


£ 

0 

P 

e,zz/Pe^0 



Unmodified 13 -moment 

Modified 13 -moment 

Exact 

1 

0.655 

4.135 

2.749 

4 

1.719 

1.719 

1.719 

OO 

1.964 

2.137 

2.513 


Also given in table VI are the modified (by information theory) values of P e z2 
as computed from equation (95) with 


V = ^ < 96) 

from equation (41). The previously mentioned possibility of overcorrections of errors 
inherent in the unmodified Grad distribution is especially apparent in the case of Coulomb 
interactions (? = 1), but evidently fails to occur when ? is greater than 4. It is inter- 
esting that Maxwellian molecules (? = 4) again have the peculiar property whereby tire 
corresponding exact distribution function of equation ( 88 ) is identical with that of equa- 
tion (93) and entropy is maximized automatically. 

Since a feature of the derivation of equation (93) from information theory is the fact 
that the square of the first-order function cp ^ contributes heavily to the resulting 
77 -dependent terms and since the < of equation (12) rather than that of equation (78) 
was used in this regard in order to preserve the 13-moment characteristics of (p Q , it 
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may be more consistent to substitute into equation (95) the following value of 77 obtained 
from the 13 -moment equations (39) and (94) rather than the “-moment value in 
equation (96): 


77 = 


5£(4 - Q 
13^ 2 - 16£ + 16 


(97) 


In particular, 


P e zz ^modified 13-moment, £= 1,77 = = 3.140p 6 /3q (98, 

represents a substantial improvement over the corresponding value in table VI, but 

P e ,zz ( modified 13-moment, £=“, 77 = - = 2.114p e 73Q (99) 


does not. 

These examples illustrate the difficulty in making general statements irrespective 
of the microscopic parameters. It does appear, however, that the type of consistency 
invoked in equations (97) to (99) is highly desirable over an important range of £ and 
is not very damaging for the remaining values of £. In any event, it seems safe to 
assert that except for values of £ in the neighborhood of 4, and especially for values 
of £ somewhat less than 4, the unmodified 13-moment approximation is quite inadequate 
for closing-out purposes at the second-order level. This trend toward greater inaccura- 
cies is expected to continue in the still higher orders, an example of which is the calcu- 
lation of the third-order contribution to the electrical conductivity. 

Reference Velocity Modification 

It was shown in the preceding paragraphs that a consistent treatment of the modi- 
fied (by information theory) 13 -moment approximation leads to a substantial improvement 
over the unmodified approximation for an important range of the interaction parameter |. 
Another modification (ref. 1) consists of taking the electron particle velocity relative to 
the mean velocity of the electrons in the Grad approximation, instead of relative to the 
mean mass velocity of the plasma. Everett's argument for this procedure is as follows: 
Although the electrons still must be reasonably close to thermodynamic equilibrium 
among themselves because the dominant term is Maxwellian relative to the electron dif- 
fusion velocity, they need no longer be close to thermodynamic equilibrium with the heavy 
particles. 
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Such an argument seems inconclusive, however, because the description of the 
frictional distortion by electron— heavy -par tide collisions of the Maxwellian distribution 
function for electrons does not appear to be generalized. Hence, the requirement of near 
thermodynamic equilibrium between species does not appear to be alleviated. The pur- 
pose of the present section is to investigate whether, in fact, Everett's modification is an 
improvement over that of information theory. 

Everett's distribution function (ref. 1) through second-order terms can be written 


f e = n e 


m, 


v 3/2 


2?rkTL 


exp 




5 

2 


( r z ' Pq) 


+ 





( 100 ) 


the primes in which indicate quantities referred to the electron frame of reference. If 
the parameters and are introduced such that 


TL 



and 


o , o 
P = P + 

^ejZz ■'e^z 


«p e Po 

15 


°2 


( 101 ) 


(102) 


the subsequent expansion of equation (100) yields 


m, 


1 3/2 


fe ° e 2irkT c 


( 1 + c o 


e, Everett J 


(103) 


and 


^e, Everett 


- 2 Mf + 1 - *») + it *•.« W - ^ + 7 


° 

2p e 


2 _ 


2( Vo L 


2 „ .2 


+ /3o( 2 Vz + «ir 2 -!«!-!) + 


4 m 


;(2r 2 - 7 ) - 


2 5 

y +0 

h 


(104) 
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The requirements that the energy 
(104) result in identities for T e and 


n 

and m e Cg Z moments of equations (103) and 
p e,zz> respectively, give 


a l = 


2 

3 


(105) 


and 


“2 = - 5 


(106) 


so that equation (104) becomes 


^e, Everett ^O^z^ 7 ^ + 1 - wj + p e,zz( 3y z 
+ ^[ y 2( 2y2 .,). y2 + |] 


(107) 


There is, indeed, a difference between 
tion of equation (12). 


^e, Everett 


and the unmodified 13 -moment func- 


If equation (107) is substituted into equation (92) for the hypothetical plasma out- 
lined previously, there results 



Everett 


5 « h §4 


+ ( 4 - S)v 

3(3ij - 2)A(^)R^gR_^ g ^ 


Pe^O 


(108) 


o 2 2 

Accordingly, the values of P g zz from equation (108) are 3.055p e /3g, 1.719p e /3Q, and 

2 9 

2.231p /3 n for £ = 1, 4, and °°, respectively, if rj is computed from equation (96). 

e u o 9 9 9 

Likewise, the values of P e zz are 2.501p e /3g, 1.719p 6 /3q, and 2.169p e /3g for 
£ = 1, 4, and °°, respectively, if 77 is computed from equation (97). 

As seen in table VII, wherein all the calculations of ?> e zz (£=l) are listed, 
Everett's function gives the best values yet of the traceless pressure tensor. It seems 
quite adequate for this purpose over the entire range of £ and is definitely to be pre- 
ferred over the distribution function derived from information theory. 

Of some concern, however, is the fact that Everett's modifications in equation (107) 
are confined to the second perturbation order; consequently, his method fails to improve 
the inaccurate entropies and relations between first-order moments that were deduced 
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TABLE VII.- TRACELESS PRESSURE TENSOR COMPONENT 
FOR COULOMB INTERACTIONS (.£ = 1) 


p e,zz/ p e% 

Method 

0.655 

Unmodified Grad 13-moment. 

4.135 

13 -moment with entropy maximization and exact heat flux. 

3.140 

13-moment with entropy maximization and approximate heat flux. 

3.055 

Everett 13 -moment with exact heat flux. 

2.501 

Everett 13 -moment with approximate heat flux. 

2.749 

Exact solution. 

previously (table 

s I to V) by using the Grad distribution. Nor can it be said without addi- 


tional calculations and comparisons that Everett's third-order results are valid. 


CONCLUSIONS 


The utility of Grad-like distribution functions has been investigated to the first and 
second perturbation orders with the following conclusions: 

1. The 13 moments in Grad's velocity distribution function are deprived of their 
character as independent variables of the problem when the function is employed in the 
Everett manner to close out the macroscopic equations of change. The relations thus 
derived between the moments are not satisfied by the exact moments because of the 
inaccurate evaluations of the collisional transfer terms using the Grad approximation. 

2. Although the errors in the Everett type of calculations of the heat flux relative 
to electron diffusion are not serious of themselves, the fact that the heat flux multiplies 
the strongest velocity polynomial in the 13-moment distribution function can lead to much 
larger errors in calculations of other quantities. 

3. The substitution into the original Grad 13-moment distribution function of the 
incorrect heat flux obtained by the closing-out method results in poor predictions of both 
the entropy and the entropy production. This is especially pronounced if the important 
difference between the collision time for entropy production and the collision time asso- 
ciated with electron diffusion is considered. That such errors are caused primarily by 
incorrect moments, rather than by the approximate form of the 13-moment function, is 
demonstrated by the excellent value of the entropy found when the correct heat flux is 
substituted into the Grad distribution. 

4. The Grad function is much more accurate in its prediction of entropy if, as orig- 
inally intended, experimental values are substituted for the 13 moments. Even so, the 
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results are quite poor under special circumstances - namely, thermodynamic and field 
conditions for which the 13 moments are insufficient to describe the frictional effects 
arising from electron — heavy -particle collisions. In no case, however, does the 
13-moment distribution function provide correct values of the moments higher than the 
heat flux. 

5. A generalization of Everett’s method to include Grad-like functions of arbitrary 
complexity shows a one-to-one correspondence between the number of moments consid- 
ered and the number of Sonine polynomials employed in the first-order Chapman-Enskog 
solution. In particular, the 5-moment analysis corresponds to the first Sonine approxi- 
mation and the 13-moment analysis corresponds to the second Sonine approximation. 
What benefit, if any, can be obtained by using the Everett method at these levels is prob- 
ably associated with the mechanics of solution. 

6. The unmodified 13 -moment distribution function of Grad is unacceptable for 

closing-out purposes at the second-order level. Two modifications were considered 
with the following result: Everett's method of referring the Grad function to the electron 

frame of reference is more successfully used in predicting the traceless pressure tensor 
than is the distribution function obtained from information theory (entropy maximization). 
Neither technique improves the first-order computations. 

7. Neither the Everett modification to the Grad function nor that of information 
theory depicts the same time dependence as the exact second-order solution derived in 
the present research. The interpretation is not altogether clear, but one must expect 
that second-order properties should depend on time in this problem through means more 
complex than the temperature. In any event, the present solution is a rather pronounced 
departure from that of Chapman and Enskog because the second -order perturbation func- 
tion appears on the left side of the second-order perturbation equation. 

8. Everett's calculations at the third and higher perturbation orders require fur- 
ther consideration before definite statements can be made regarding their validity. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., June 23, 1969. 
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APPENDIX A 


DEVELOPMENTS FROM INFORMATION THEORY 


It is shown in this appendix that the first-order Grad distribution function of equa- 
tion (14) can be derived just as well from the concepts of information theory as from the 
original Grad expansion in terms of three-dimensional Hermite polynomials (ref. 1). A 
start in this direction was recently achieved by Stankiewicz (refs. 4 and 5), but he stopped 
short of the required number of moments. 

The fundamental postulate of information theory states that the best distribution 
function that can be derived from limited experimental knowledge alone is the one which 
maximizes the entropy density 



subject to the measured constraints, 
consist only of measurements of the 
method of Lagrange multipliers jl j 
tional expressions: 


If the nonequilibrium contributions to the latter 
£j defined in equation (3), the application of the 
to 6s = 0 and 6^ = 0 yields the following varia- 


and 



(A2) 


<Pe 


I 


y 



* y 


(A3) 


The combination of equations (3) and (A3) yields 


\ | £ 2_j " k (2j + 2k + 3): : j a. (A4) 

j 

the substitution of which into equation (A3) gives the 16-moment perturbation function 

|p - f)(?i - 1 h) ■ y * k( r * - ^ * “)(% - 7 % + f ? o) • ? (as) 

if j assumes the values 0, 1, and 2. The description of this function as a 16-moment 
approximation refers to the addition of ^ to the original 13 moments of Grad. 
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APPENDIX A - Concluded 


Among other requirements of a similar nature, the condition 

^2 = 7 ^1 0 O ) ( A6 ) 

must be satisfied before Grad's 13-moment function is consistent to first order with 
information theory and the full complement of ^ constraints. 
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APPENDIX B 


COLLISION INTEGRALS 


It is shown in this appendix that equation (83) follows directly from the collision 
dynamics when the electron — heavy -particle interaction potential is written as the 
inverse £-power of the separation distance. 

The first step in such a development is the formulation of the reduced electron par- 
ticle velocity after a collision y' in terms of the reduced electron particle velocity 
before the collision y and the corresponding scattering or deflection angle y. This 
straightforward exercise in geometry yields 

y' = y cos x + 7 sin y[[ (sin cos e + cos 9 cos sin e) 

- j(cos cp cos e - cos 9 sin c p sin e) - k sin 6 sin el (Bl) 


where 9 and cp are the polar and azimuthal angles, respectively, of y in spherical 
coordinates and e is the azimuthal angle of y’ with respect to the direction of y. 

Accordingly, 



b db de 


2rry \ (1 - cos y)b db 

J 0 


(B2) 


The remaining integral involves, of course, the detailed collision dynamics, but that cal- 
culation is circumvented here by the following use (ref. 2) of equations (15) and (35): 
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(B3) 


and 

f - b db de = 2,f (1 - cos X )b db = 7 , -VI- ( B4) 

where nj is the number density of the j-type heavy particles. 
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APPENDIX B - Concluded 


Finally, 


j* (yl - r; 2 ) b db de = J y 2 (l - cos 2 y) + 2 y L 


'sm x cos X 


is 



= it ^ 3 y 2 ~ y2 )J 0 ( J _ cos2 )c) b db 
= A(£)7^3y 2 - y 2 ) J (1 - cos y)b db 
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(B5) 


from equations (Bl) and (B4). The ratio A(£) of the two impact parameter integrals 
has the values 2, 1.034, and 0.666, respectively, for £ = 1, 4, and <*>. (See ref. 6, 
pp. 546-549.) 


35 



REFERENCES 


1. Everett, Willis L.: Generalized Magnetohydrodynamic Equations for Nonequilibrium 

Plasma Systems. Rarefied Gas Dynamics, Vol. I, J. A. Laurmann, ed., Academic 
Press, 1963, pp. 1-25. 

2. Meador, Willard E.: A Semiempirical Collision Model for Plasmas. NASA TR R-310, 

1969. 

3. Chapman, Sydney; and Cowling, T. G.: The Mathematical Theory of Non-Uniform 

Gases. Second ed., Cambridge Univ. Press, 1952. (Reprinted 1960.) 

4. Stankiewicz, Norbert: Thermodynamic Study of Plasmas Using the Principle of 

Maximum Entropy. NASA TN D-4621. 1968. 

5. Stankiewicz, N.: The Ellipsoidal Velocity Distribution Function. NASA TN D-5026, 

1969. 

6. Hirschfelder, Joseph O.; Curtiss, Charles F.; and Bird, R. Byron: Molecular Theory 

of Gases and Liquids. John Wiley & Sons, Inc., c.1954. (Reprinted with corrections 
1964.) 


36 


NASA -Langley, 1969 25 L-6523 



